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Abstract 



A theory is presented for the frequency dependence of the power spectrum of 
photon current fluctuations originating from a disordered medium. Both the 
cases of an absorbing medium ("grey body") and of an amplifying medium 
("random laser") are considered in a waveguide geometry. The semiclassi- 
cal approach (based on a Boltzmann-Langevin equation) is shown to be in 
complete agreement with a fully quantum mechanical theory, provided that 
the effects of wave localization can be neglected. The width of the peak in 
the power spectrum around zero frequency is much smaller than the inverse 
coherence time, characteristic for black-body radiation. Simple expressions 
for the shape of this peak are obtained, in the absorbing case, for waveguide 
lengths large compared to the absorption length, and, in the amplifying case, 
close to the laser threshold. 

PACS: 42.50.Ar, 05.40.-a, 42.68.Ay, 78.45. +h 



I. INTRODUCTION 

The noise power spectrum of a black body is frequency independent for frequencies below 
the absorption band width. The inverse of the band width is the coherence time r co h of the 
radiation 10, which for a black body is the longest relevant time scale — hence the white 
noise spectrum P(fl) for Q < l/T coh . In a weakly absorbing, strongly scattering medium 
there appear two longer time scales: The absorption time r a and the time L 2 / D it takes to 
diffuse (with diffusion constant D) through the medium (of length L). As a consequence, 
-P(fi) for such a weakly-absorbing medium (sometimes called a "grey body") starts to decay 
at much lower frequencies than for a black body having the same coherence time. 

Although there is by now a substantial literature on the theory of grey-body radiation 
the results have been limited to either the zero or high-frequency limits of the noise 
spectrum (or, equivalently, to short or long photodetection times). In the present work 
we remove this limitation, by computing P(fl) for a diffusive medium for arbitrary ratios 
of Q, 1/T a , and D/L 2 . We compare two different approaches in a waveguide geometry: 
One which is fully quantum mechanical (based on random- matrix theory |7]||) and another 
which is semiclassical (based on a Boltzmann-Langevin equation |]]). Each method has 
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its advantages and disadvantages: The quantum theory includes interference effects, which 
are ignored in the semiclassical theory, but it is mathematically more involved. Complete 
agreement between the two approaches is obtained in the limit that the waveguide length L 
is much smaller than the localization length (equal to the mean free path times the number 
of propagating modes). 

The results for absorbing media can be applied directly to linear amplifiers, by formally 
changing the sign of the temperature and the absorption time. Loudon and coworkers pH|JTT|] 
used this relationship to calculate the noise power spectrum of a waveguide without disorder. 
The generalization to a diffusive medium presented here describes a random laser [|T2|] below 
threshold. 

The outline of this paper is as follows. We start with the semiclassical approach, pre- 
senting a general solution of the Boltzmann-Langevin equation in Sec. II and applying it to 
a waveguide geometry in Sec. III. The quantum mechanical approach is developed in Sec. 
IV. For the quantum theory we need the correlator of reflection and transmission matrices 
at different frequencies. These are calculated in the appendix, using the random-matrix 
method of Ref. |13j. We discuss our findings in Sec. V. 



II. SEMICLASSICAL THEORY 



Starting point of the semiclassical theory is the Boltzmann-Langevin equation for photons 
of Ref. |§. We first consider an absorbing medium (in equilibrium at temperature T), leaving 
the amplifying case for the end of this section. We make the diffusion approximation, valid 
if the mean free path I is the shortest length scale in the system (but still large compared 
to the wavelength). The fluctuating number density n(uj,r, t) and current density j(u;,r, t) 
of photons at frequency lu, position r, and time t are related by |j| 

i = -D^ + C u (2.1) 
or 

- ¥t + Wv -i = DC 2 { P f-n) + C«. (2.2) 

Here D = |cZ is the diffusion constant, £ a = \JDr a is the absorption length (with r a the 
absorption time), p = Amuj 2 (2ttc)~ 3 is the density of states (not counting polarizations), and 
/ = [exp (Hu/kT) — is the Bose-Einstein function. We assume £ a ^> I. The fluctuating 
source terms Cq and C\ have zero mean and correlators 



Cq(w, r, t)£ Q (u', r', f) = 5(lu - u')6(t - t')5(r - r')D^ 2 (2fn + pf + n), (2.3a) 
C la (LU,r,t)C 1/3 (uj',v',t') = 25 a(3 5(tu - u')5(t - t')6{r - v')Dn{l + n/p). (2.3b) 

The cross-correlator of £ an d C>\ is given in Ref. 0, but will not be needed. Combining 
Eqs. (|2.1|) and ( p.2|) we find equations for the mean n and the fluctuations Sn of the photon 
number density n = n + 5n, 

1 dn d 2 n n pf 

1 dSn d 2 Sn dn 1 d Co 

+ ^9 ~F> = TT^T ' £i FT- ( 2 -5) 



D dt dr 2 e. Ddr 1 D 
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FIG. 1. Thermal radiation (solid arrow) is incident through port So on an absorbing disordered 
medium (shaded). The outgoing radiation (dashed arrows) is absorbed by photodetectors. 



We present a general solution for the multiport geometry of Fig. 1. Thermal radiation 
is incident through the port S and can leave the system via ports S , Si, S 2 , ■ ■ ., where it 
is absorbed by photodetectors. The corresponding boundary conditions are n(u, r, t)\ re s p — 
n in (uj,t)Sp . We assume that the closed boundaries £ of the system (with volume V) are 
perfectly reflecting. The separation of the ports is of order L ^> I. In what follows we assume 
detection of outgoing radiation in a narrow frequency interval 5u around to. We require that 
Sou is small both compared to uj and to l/r^h. To minimize the notations in this section we 
omit the frequency argument uj and use units in which 8u = 1. (We will reinsert Soj in the 
next section.) 

The Green function of the differential equations ( |2.4| ) and ( |2.5| ) in the Fourier represen- 
tation with respect to the time argument satisfies 

(^-C 2 + §)G(r,r',fi) = 5(r-r'). (2.6) 

[Fourier transforms are denned as /(fi) = /f^ dt e int f(t).] For frequency resolved detection 
we require Q 5uj. We impose the boundary conditions 

G(r,r',fi)| rG5p = 0, p = 0,1,2,..., (2.7a) 

s . a^s. Ue = o, ( 2 .7b) 

or 

where S denotes the outward normal direction to the surface S. We consider separately the 
mean and the fluctuations of the photon number and current densities. 

A. Mean solution 

The average photon density satisfying Eq. fl2.4p can be expressed in Fourier representation 
in terms of the Green function fl2.6|) . 
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n(r, ft) = -2npfC 2 m J & G(r, r', 0) + n in (ft) | dS' • ^j', Q) . (2.8) 

V So 

Substituting this formula into the expression for the current ( |2.1| ) and integrating over the 
area S p one obtains the mean outgoing current I p through port p ^ 0, 



J p (ft) = 2tt Pj D/C 2 <H^) JdS-Jdr 



, dG( T,T',0) 

dr 

, O z G(r,r'Sl) 
i 



v 

}2 



ZM in (ft) J dS a J dS' / . (2.9) 

5 P 5o 



(Summation over the repeating Greek indices is implied.) The first term oc <5(ft) is the 
time-independent mean thermal radiation from the medium. The second term is that part 
of the mean radiation entering through port that leaves the medium through one of the 
other ports. (The restriction to p ^ is not essential but simplifies the general formulas 
considerably, so we will make this restriction in what follows.) 



B. Fluctuations 



The fluctuations in the number density follow in a similar way from the Green function 
and Eq. (Pfl, 

5n(v, ft) = i J dr' G(v, r', ft) ( A • C^v', ft) - £ (r', ft)) + 5n in (ft) J dS' • ^ . 

(2.10) 

The fluctuation of the current density is then given by Eq. ( |2.1| ), 

6j a {r,Q) = J dr' (G a p(r,r',{l)£ w (r',n) + ^ C (r', ft)) 

- D5n in (Q) J dS'p G Q/3 (r,r',ft). (2.11) 

So 

We have defined 

Gae(r,r',Q) = 92 f^ ] +V(r-r'). (2.12) 
We seek the correlator of the current fluctuations 



C Q/3 (r, ft; r', ft') = 6j a (r, ft)^(r', ft) (2.13) 



for r G 5 p , r' G with p, g ^ 0. With the help of Eqs. fl2.3|) and ( |2.11| ) it can be expressed 
as 
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C Q/3 (r,fi;r,fi) = — dr — [(2/ + l)n(r ,Q + Q)+pf] 

Sa £ ° r a ar j3 

r r 1 /■ r/Q" 
+2D / dr"G Q7 (r, r", fi)GWr', r", fi') n(r", Q + Q') + - n(r", + fi")n(r", Q' - SI") 

J l O J 2lX 

v r 

(2.14) 

Following Ref. |J, we have neglected the term oc 8n in in Eq. ( |2.11| ) (smaller by a factor l/L) 
and the cross-correlator C$Ci (smaller by a factor Z/£ a ). 

We now integrate r and r' over S p and S 9 to obtain the correlator of the total currents 
through ports p and q, 

c pq (n,Q') = J ds a J ds' p c aP {v,n-v',n') = <7#(n,tf) + cg>(n,n'). t 2 -^) 



The first term contains the contribution from the terms linear in the number density 



n 



in Eq. (|2.14 ). Performing integration by parts and using Eqs. ( p . 6| ) — ( p78|) we find that this 
term vanishes for p ^ q. For p — q it contains the mean current, 

cW(n,n') = 6 pq i p {n + n'). (2.16) 

For a time-independent mean current I p one has a white-noise spectrum (7^(0,°/) = 
27T(5pq5(fi + Q!)I P . This is the usual shot noise, corresponding to Poissonian statistics of 
the current fluctuations. The second term C-® describes the deviations from Poissonian 
statistics. It arises from terms in Eq. (|2 . 1 4j) that are quadratic in n. Performing again an 
integration by parts, one finds 

Sp Sq V 

x dG(r,T",n) dG(v',v",n') 

dr a dr'p 

Equation (|2.17| ) together with Eq. (|2.8|) is the result that we need for our analysis of the 
frequency dependence of the noise spectrum. 



C. Amplifying medium 

The extension of our general formulas to an amplifying medium (in the linear regime 
below the laser threshold) is straightforward 0: We assume that the frequency uj at which we 
are detecting the radiation is close to the frequency of an atomic transition with (on average) 
^upper and iV lower atoms in the upper and lower state. Then the Bose-Einstein function can 
be replaced by the population inversion factor / = N uppeT (Ni oweT — A^ upper ) _1 . This factor 
is negative in the amplifying case (when iV upper > A^ ower ), with / = — 1 for a complete 
population inversion. ( Equivalent ly, one can evaluate / at a negative temperature |lTJ , with 



T — > 0~ for complete inversion.) An amplifying medium has a negative absorption time 
r a = il/D. We can account for this by taking £ a imaginary. With these two substitutions 
for / and £ a our formulas for an absorbing medium carry over to the amplifying case. 
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FIG. 2. Thermal radiation (solid arrows) is incident on a waveguide containing an absorbing 
or amplifying disordered medium. The transmitted radiation (dashed arrows) is absorbed by a 
photodetector. 

III. WAVEGUIDE GEOMETRY 

For the application of our general formulas we consider a waveguide geometry (see Fig. 
2). The waveguide has length L and cross-sectional area A, corresponding to N = u 2 A/47rc 2 
propagating modes (not counting polarizations) at frequency uj. We abbreviate s = L/£ a . 
We consider a stationary incident current J = ^cA5un in = (N5u/27ip)n in , and calculate 
the noise power spectrum of the transmitted current, 

oo 

P{Q) = J dt e int 6I{t)6I(0). (3.1) 

— oo 

In terms of the correlator of the previous section, one has Cu(Q, Q') = 2nP(Q)5(fl + Q'). 



A. Absorbing medium 

We calculate the noise power from Eqs. ( |2.8|) and ([2.17|) , using the Green function 



sinh [(x < /^ a )y/l - iVtT a \ sinh [(s - x > /^ a )y/l - iVtr a 



sinh [si/1 — i£lT a ] 



(3.2) 



where x< and x> are the smallest and largest of x, x', respectively. The mean photon density 
is time independent. In Fourier representation one has, from Eq. ( [2.8|) , 



n\x. 



, Vl) = 2n5{Vl) — - — (sinh s — sinh (x/t; a ) — sinh (s — x/£ a 
sinh s v 

, o xfn\- sinh(s-x/^ a ) 

+ 27tb(il)n- m — . 

sinh s 



(3.3) 



The mean current I = I t h + /trans is the sum of the thermal radiation from the medium 



ADf 

J th = ^-(iV(W27r)tanh(s/2 



and the transmitted incident current 

I 



trans 



c£ a sinh s 
6 



(3.4) 



(3.5) 



Substitution of Eqs. ( |3.2|) and ( |3.3|) into Eq. (|2.17p yields the super-Poissonian noise 
P — I as a sum of three terms, P — I = P th + Pt rans + Pex> with 

^ 8Pf 2 Tr . . ? . /cosh (s — s') — coshs'\ 2 , w . . 

Pth(fi) = -^-{N8u/2<k) J ds' y 1 ginh ; g J K(s',s), (3.6) 

SDf 2 f cosh 2 fs - s"l 

P tra „s(fi) = ^(2n/N5u) / ds' | g 2 Sj K(^, S ), (3.7) 

c^a ^ sinn s 

Pex(^) = /^ [COShg/ ~ COSh(g " 2 S/)]COSh(g " g0 ^ g ). (3.8) 



We have defined 



K(s', 



sinh ( s' \/l — iVtT a ) 



(3.9) 



sinh ( s a/1 — iVtT a ) 

The two terms Ptrans and P tn describe separately the noise power of the transmitted incident 
current and of the thermal current from the medium. The term P ex is the excess noise due 
to the beating of the incident radiation with the thermal fluctuations from the medium. 

The three contributions are plotted separately in Fig. 3. For L ^> £ a the frequency 
dependence simplifies to 

Pth(n) = (3.io) 

1 + x 

Prans(fi) = ^^(2tt/N5u) + ^ > C 3 ' 11 ) 

loiJ \ x — 1 X z + X J 

1 + 2x 

Pex(^) = //trans ; 95 (3.12) 

X + X 2 

where we have defined 



x = Re yl - zfir a = [|(1 + ^V a 2 ) 1/2 + |] 1/2 . (3.13) 

As discussed in Ref. || (for the zero-frequency case) the result for Ptrans requires that the 
incident radiation is in a thermal state, at some temperature To. (The quantity f(u,To) = 
Iq{2tt / N8uj) is the corresponding value of the Bose-Einstein function.) There is no such 
requirement for P tn and P ex , which are independent of the incident state. For T > T we 
may generally neglect P tn and P ex relative to P tr ans, so that P = /trans + Ptrans- However, if 
the incident radiation is in a coherent state, then Ptrans = and since for sufficiently large Iq 
we may neglect P t h, we have in this case P = /trans + P ex - The contribution P tn is important 
mainly in the absence of external illumination, when P = / tn + P tn . 



B. Amplifying medium 

The results for an amplifying medium are obtained by the substitution £ a — > i£ a , f — > 
A upP er(A lower - Auppe,.)- 1 , cf. Sec. IIC. The frequency dependence of P th , Ptrans, and P cx 
following from Eqs. (|3.6|) - (|3.8|) is plotted in Fig. 4 for lengths L below the laser threshold 
at L = 7r£ a . 
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FIG. 3. Frequency dependence of the three super-Poissonian contributions to the noise power, 
P — I = Pth + -ftrans + Pex, for different values of s = L/^ a in an absorbing waveguide. The 
solid curves are computed from Eqs. (|3.6|)-( j3~8l ), the dashed curves are the large- s asymptotes 
(|3l0|)-(|3~ll). The parameter Z is defined as Z = (c(, a /2D)(2ir/N5uj). 



C. Cross-correlator 



In the absence of any incident radiation, the noise P — Jth + -Pth is due entirely to 
the thermal fluctuations in the medium. The current fluctuations at the two ends of the 
waveguide are correlated, as measured by the cross-correlator 



Pia(fi) = / dt e int 5h{t)5I 2 (0). 



(3.14) 



From Eqs. (|2~T7| ), 0, and Q we obtain 

_ 8Df 2 .^ TC /cosh (s - s') - cosh s r 

P 12 (s i) = _^ ( i Vfo /2.)/ rfs ' ( sinh ' s 
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FIG. 4. Same as Fig. 3, for the case of an amplifying waveguide. The laser threshold occurs at 

S = TT. 



sinh[s'\/l ~~ i^T a ] sinh[(s — s')\/l + iQr a 



(3.15) 



The cross-correlator is plotted in Fig. 5 for both the absorbing and amplifying cases. The 
outgoing currents at the two ends of the waveguide are anti-correlated for flr a ;§> 1. 



IV. COMPARISON WITH QUANTUM THEORY 

A fully quantum mechanical theory for the photocount distribution of a disordered 
medium was developed in Refs. In this section we verify that it agrees with the 

semiclassical results of the previous section. We consider the same system of Fig. 2, a dis- 
ordered waveguide with a photodetector at one end and a stationary current incident at the 
other end. We assume that the incident current originates from a thermal source at tem- 
perature Tq. The photocount distribution is the distribution of the number of photons n(t) 
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FIG. 5. Frequency dependence of the cross-correlator of the outgoing current at the two ends 



of the waveguide, in the absence of any external illumination. Computed from Eq. (3.15) for the 
absorbing case (lower panel) and amplifying case (upper panel). 

counted (with unit quantum efficiency) in the time interval (0, t). Substitution of I — dn/dt 
in the definition ( pM| ) of the noise power P(£l) leads to a relation with the variance Varn(t) 
of the photocount, 



P(Q) = -VI 2 J dtVam(t) cos fit, 
o 

Var nit) = [ dQ Q- 2 P{Q) (cos Qt - 1) . 

7T J 



(4.1a) 
(4.1b) 



The variance can be separated into two terms, Varn(t) = n(t) + n{t) = tl + K,(t), with 
K{t) the second factorial cumulant. The term tl, substituted into Eq. (|4.1a|) , gives the 
frequency-independent shot noise contribution I to the power spectrum, 



P(Q) =1 -Q 2 J dt K(t) cos Qt. 
o 



(4.2) 



The cumulant k = K trans + Ktt + K ex contains separate contributions from the transmitted 
incident radiation and thermal fluctuations in the medium, plus an excess contribution from 
the beating of the two. These contributions have an exact representation in terms of the 
N x N reflection and transmission matrices r(u), t(u) of the waveguide MM, 
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oo , oo 



^ans(t) = J^J ^L{u-J,t)f{u],T Q )f{u',T Q )TrT{u)T{u'), (4.3) 



= hr ~ *) T )^'' T ) Tr <2HQ(o/), ( 4 - 4 ) 

J 2ix Jo 2tt 
o 

^ex(t) = / ^ - 2/(0;, T )/(a/, T)Tr T(co)Q(u;'), (4.5) 





where we have defined 

t 



L(u,t) = Jdt' J dt" exp[iu;(t' - t")} = 2u~ 2 (l - coswt), (4.6) 



Q{u) = II - r(uj)r\u) - fH^(w), (4.7) 

T(uj) =t(uj)t\uj). (4.8) 



Substitution into Eq. (|4.2j) gives the corresponding contributions to the noise power 

P — I -\- -Ptrans + -Ph + -^ex, 

duj 



Ptrans(fi) = ^ / ^/(w, T )f{u + ft, T )Tr T( W )T( U + ft) + {ft - -ft}, (4.9) 



P th (ft) = -j JHf{u, T)f(u + ft, T)Tr Q{uj)Q{uj + ft) + {ft - -ft}, (4.10) 

2 J ZTT 


P cx (ft) = - / -^2/(o;, T )f{u + ft, T)Tr r(w)Q(o; + ft) + {ft -+ -ft}. (4.11) 

ZJ Z7T 


As in the previous section, we assume a frequency-resolved measurement in an interval 
5uj <C ou, l/r co h with ft <C 5u. We may then omit the integral over uj and approximate the 
argument uj ± ft in the functions / by a;. We take the ensemble average (■ • •) of the noise 
power, in which case the contributions from ±ft are the same. Finally, we insert the incident 
current Jo = f(u>,To)N5uj/27r, to arrive at 

Ptrans(ft) = (2n/N5uj)I 2 (N- 1 TrT(uj)T(uj + ft)), (4.12) 
P th (ft) = (^/2vr)/ 2 (^,T)(Ar- 1 Trg(cu)g(cu + ft)), (4.13) 
P cx (ft) = 2/ /(^,T)(iV- 1 TrTHg(a; + ft)). (4.14) 

It remains to evaluate the ensemble averages. This is done in the appendix, by extending 
the approach of Ref. |l3j to correlators of reflection and transmission matrices at different 
frequencies. The calculation applies to the diffusive regime that the length L of the waveguide 
is large compared to the mean free path I, but still small compared to the localization length 
Nl. (The absorption length £ a is also assumed to be 3> I.) The results are 



11 



(Ar- 1 Trr(cu)T(c; + fi)) 



8D 



8D 



ds' K(s', s 



cosh (s — s') 



(N- lr Tr Q{u)Q{u + tt)) = — I ds' K(s', s 



sinh s 

, [cosh s' — cosh(s — s')} 2 



sinh s 



(N- lr TrT(oj)Q(uj + ty) = — / ds'K(s',s) 



, . cosh(s — s') cosh s' — cosh 2 (s — s') 



sinh s 



(4.15) 
(4.16) 
(4.17) 



where s = L/^ a and the kernel K(s', s) is defined in Eq. (|A29 ). The combination of Eqs. 
( |4.12j )- (|4.17| ) agrees precisely with the results ( p.6|) -( ^8|) of the semiclassical theory. The 
quantum theory is more general than the semiclassical theory, because it can describe the 
effects of wave localization. The method of Ref. [13| gives corrections to the above results 
in a power series in L/Nl. We will not pursue this investigation here. 



V. DISCUSSION 



We have presented a theory for the frequency dependence of the noise power spectrum 
.P(f2) in an absorbing or amplifying disordered waveguide. The frequency dependence is 
governed by two time scales, the absorption or amplification time r a and the diffusion time 
L 2 /D, both of which are assumed to be much greater than the coherence time r co h. A 
simplified description is obtained, in the absorbing case, for lengths L much greater than 
the absorption length £ a = \JDr a , and, in the amplifying case, close to the laser threshold 
at L — 7r£ a . We will discuss these two cases separately. 



A. Absorbing medium 



The general formulas (P7B|)— (P7B|) for P = I + P th + P trSkns + P cx simplify for L > £ a to 
Eqs. (|3.10|) - (|3.12|) . To characterize the frequency dependence we define the characteristic 



frequency VL C as the frequency at which the super-Poissonian noise has dropped by a factor 
of two: 

p(n c )-J=J(p(o)-j). (5.i) 

In the absence of any external illumination (J = 0) we have, from Eq. ( p.l(J| ), 

P = /th(l + -f— J, /th = ^r-{N8u /2tt), (5.2) 



with x = Re a/1 — iVLr a} hence Vt c = 17 /r a . If the illumination is in the coherent state from 
a laser, then we have, from Eq. ( |3.12| ), 

P = Arans ( 1 + / 2 I , Arans = — T~ e , (5.3) 
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FIG. 6. Ratio of / 2 n and -Pth in an amplifying waveguide as a function of its length for different 



frequencies, computed from Eqs. (3.4) and (3.6). The approximation ( |5.5[) valid near threshold for 
small frequencies is shown dashed. 



here Q c = 9/r a . In both these cases the diffusion time does not enter in the frequency 
dependence. This is different for illumination by a thermal source at temperature To much 
greater than the temperature of the medium. From Eq. ( |3.11| ), with f = f(u,T Q ), we then 
have 



P trans (^) Arans ( 1 ~t~ ^ ^ 



1 _ e -2*(x-l) 3x + 2 
2T + X 



X 



1 



(5.4) 



The characteristic frequency f2 c = (64D/L 2 r^) 1 / 4 now contains both the diffusion time and 
the absorption time. 



B. Amplifying medium 

In the amplifying case the noise power becomes more and more strongly peaked near zero 
frequency with increasing amplification. Close to the laser threshold at s = ir the frequency 
dependence of P tn for small frequencies Qr a <C 1 has the form 

p = ^th j = A f (r r\ 

th 27r[ft 2 T a 2 + 4(l-s/7r) 2 ]' th Z(tt-s)" 1 J 

Here again Z = (c^ a /2D)(2-n-/N5uj). Close to threshold the peak in the noise power spectrum 
has a Lorentzian lineshape with half-width Q c = (2/r a )(l — L/7r£ a ). At the laser threshold 
both P th and 7 t h diverge, but the ratio 7 2 h /P t h remains finite (see Fig. 6). 

Finally, we note the fundamental difference between the time scales appearing in the noise 
spectrum for photons, on the one hand, and electrons, on the other hand. The absorption or 
amplification time r a obviously has no electronic analogue. The diffusion time L 2 /D appears 
in both contexts, however, the electronic noise spectrum remains frequency independent for 
Q > D/L 2 ||14|| . The reason for the difference is screening of electronic charge. As a result 
the characteristic frequency scale for electronic current fluctuations is the inverse scattering 
time D/l 2 , which is much greater than the inverse diffusion time D/L 2 . 
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APPENDIX A: CORRELATORS OF REFLECTION AND TRANSMISSION 

MATRICES 

To compute the noise power spectrum in the quantum mechanical approach of Sec. V, we 
need the correlators of reflection and transmission matrices t(u±) and r(cu±) at two different 
frequencies u± = uo ± Q/2. (For O < w this is the same as the correlator at frequencies to 
and u + Q.) We calculate these correlators for a waveguide geometry in the diffusive regime, 



by extending the equal- frequency (Q = 0) theory of Brouwer [13 . 

Upon attachment of a short segment of length 8L to one end of the waveguide of length 
L, the transmission and reflection matrices change according to 

t ^ t 5L (l + rr 5L )t, (Ala) 
r^r' SL + t SL {l + rr SL )rt] L , (Alb) 

where the superscript T indicates the transpose of a matrix. (Because of reciprocity the 
transmission matrix from left to right equals the transpose of the transmission matrix from 
right to left.) The transmission matrix t§L of the short segment at frequency u± may be 
chosen proportional to the unit matrix, 

/ 8L 8L itt8L\ IK . 

t SL = 1 - — ± 1. A2 

V 2V 2c'r a 2d J K J 

The mean free path I' = 41/3 and the velocity d = c/2 represent a weighted average over 
the N transverse modes in the waveguide. 

Unitarity of the scattering matrix dictates that the reflection matrix from the left of 
the short segment is related to the reflection matrix from the right by r' SL = —r\ L . We 
abbreviate = Sr. The matrix Sr is symmetric (because of reciprocity), with zero mean 
and variance 

(MO — (N + l)-\5 km 5 ln + 8 kn 8 lm )8L/l'. (A3) 
The resulting change in the matrix products ttf and rr^ is 

tt f -> (1 - SL/i - 5L/c'T a )tt^ + {r5rt){r5rty 

+ r5rtt ] + (r5rtt ] )\ (A4a) 
rr^ — > (1 — 25L/1' — 25L/cT a )rr^ + (r8rr)(r8rr)^ + 8r'8r 

+ rSrrr^ + (rSrrr^)^ — r8r — (r8r). (A4b) 
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The frequency f2 does not appear explicitly in these increments. 
We define the following ensemble averages 



-R, = (JV _1 Tr (1 -rr f )), (A5) 
C = (iV~ 1 Tr (1 — r_r+)}, (A6) 
T=(N- 1 Trttf), (A7) 

where r, t are evaluated at frequency to and r±, t± at frequency a>±f2/2. Similarly, we define 
the correlators 

C rr = (jV -1 Tr (11 - r_rl)(l - r + r+_)), (A8) 

C rt = (N^Tt (11 - r_r^)t + 4), (A9) 

C tt = (N- lr Irt_tlt + t\.). (A10) 

We will see that, in the diffusive regime, these 6 quantities satisfy a coupled set of ordinary 
differential equations in L. 

The diffusive regime corresponds to the large- N limit, in which the length L of the 
waveguide is much less than the localization length Nl. In this limit we may replace Eq. 
( |A~3| ) by (Sr^Sr^) = (5L/Nl')8k m Si n . In the large-A?" limit we may also replace averages 
of products of traces by products of averages of traces. From Eq. (|A4j) we thus obtain the 
diffential equations 

l' d ^- = 2 1 {l-TZ)-TZ\ (All) 

l 'jL = 27(1 + ^ Ta)(1 ~ C) ~ c2 ' (A12) 

l 'lk = ~ lT ~ (A13) 
dC 

t ,„ = +c + c * + 2K)C rr + 2TZ{n + 2 7 ), (A14) 
dL 
dC 

I'—jy- = -(3 7 + C + C*+ K)C rt - TC rr + 2(71 + 7 )T, (A15) 

« = _(2 7 + c + C*)C tt - 2TC rt + 2T 2 , (A16) 
dL 

with the definition 7 = l'/c'r a . The initial conditions are that each of these 6 quantities — > 1 
for L -> 0. 

This set of differential equations may be simplified further if we assume, as we did in 
the semiclassical theory, that the mean free path is small compared to both the absorption 
length and the length of the waveguide. All 6 quantities ( |A5f) - (|Al0|) are of order y/j, which 
is <C 1 if I' <C c'r a , so that we obtain in leading order 

^ = 2 7 -^ 2 , (A17) 

/'^ = 2 7 (l+^r a )-C 2 , (A18) 
dL 
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^ dCf 



i 



dL 

■fdC r t 

~dL 

,dCtt 



dL 



-(C + C* + 2K)C rr + 2ft 2 , 

-{C + C* + TZ)C rt - TC rr + 2TZT, 

-(C + C*)C tt -2TC rt + 2T 2 . 



(A19) 
(A20) 
(A21) 
(A22) 



As initial condition we should now take that the product of each quantity with L remains 
finite when L — > 0. 

Although the differential equations are coupled, they may be solved separately for 1Z, C, 
T, C rr , C rt , C tt , in that order. In terms of the rescaled length s = {2 r )) l / 2 L/V = L/£ , the 
results are 



K 

C 

T 



(27) 



1/2 



tanhs 

(2 7 ) 1 /yi + ^r a 
tanh sy/1 + iVtr a 

(2 7 ) 1/2 ; 
sinhs ' 

3 7 ) 1/2 " 



C rr = a — / " s I s j s ) cosh s , 

sinh s J 





1/2 s 



(87) 1 f 

C r t = 2 — / ds' K{s , s) cosh(s — s ) coshs', 

sinh s •/ 


C tt = / ds' K{s', s) cosh 2 (s - 8 % 

sinh s 7 



(A23) 
(A24) 

(A25) 
(A26) 

(A27) 

(A28) 



where the kernel K is defined by 



K{s',s) 



sinh s'yl + i£lr a 



sinh syl + iVtr a 



(A29) 



These are the expressions used in Sec. 4 (where we have also substituted 1/27 = 4_D/c£ a ). 
The remaining integrals over s' may be done analytically, but the resulting expressions are 
rather lengthy so we do not record them here. For Q = our results reduce to those of 
Brouwer |13j (up to a misprint in Eq. (13c) of that paper, where the plus and minus signs 
in the expression between brackets should be interchanged). 
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